Intermittent magnetic field excitation by a turbulent flow of liquid sodium 
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The magnetic field measured in tlie Madison Dynamo Experiment shows intermittent periods of 
growth when an axial magnetic field is applied. The geometry of the intermittent field is consistent 
with the fastest growing magnetic eigenmode predicted by kinematic dynamo theory using a laminar 
model of the mean flow. Though the eigenmodes of the mean flow are decaying, it is postulated 
that turbulent fluctuations of the velocity field change the flow geometry such that the eigenmode 
growth rate is temporarily positive. Therefore, it is expected that a characteristic of the onset of a 
turbulent dynamo is magnetic intermittency. 
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Determining the onset conditions for magnetic field 
growth in magnetohydrodynamics is fundamental to un- 
derstanding how astrophysical dynamos such as the 
Earth, the Sun, and the galaxy self-generate magnetic 
fields. These onset conditions are now being studied in 
laboratory experiments using flows of liquid sodium [1] . 
The conditions required for generating a dynamo can be 
determined by solving the magnetic induction equation 

— =i?mV X V X B + V^B, (1) 
ot 

where B is the magnetic field, V is the velocity field 
scaled by a characteristic speed Vq, and the time is scaled 
to the resistive diffusion time = /xqctL^. The mag- 
netic Reynolds number is Rm = fioaLVo, where a is 
the conductivity of the fluid and L is a characteristic 
scale length [2]. In the kinematic approximation, the 
velocity field is assumed to be a prescribed flow (either 
the flow is stationary or its time dependence is speci- 
fied) and Lorentz forces due to the magnetic field are 
neglected. Equation 1 is then linear in B and can be 
solved as an eigenvalue problem. Several different types 
of stationary, helical flows have been shown theoretically 
to be kinematic dynamos [3, 4, 5, 6], which in turn has 
lead to the design of current dynamo experiments. For 
particular flows, the kinematic model predicts a critical 
value of the magnetic Reynolds number, Rrricrit, above 
which the magnetic field becomes linearly unstable, i.e. 
for Rm > Rrricrit a small seed magnetic field will grow ex- 
ponentially in time. The dynamo onset conditions have 
been tested in helical pipe-flow experiments at Riga [7, 8] 
and Karlsruhe [9, 10]. Both experiments generated mag- 
netic fields at a value of Rrricrit consistent with predic- 
tions from the kinematic theory. 

Fluids and plasmas such as the Earth's liquid core, 
the solar convection zone, and liquid metals are turbu- 
lent under the conditions required for a dynamo. The 




FIG. 1: A schematic of the Madison Dynamo Experiment 
with superimposed magnetic field lines of the theoretically 
predicted dominant magnetic eigenmode. 

Riga and Karlsruhe experiments use highly constrained 
flows to create the helical geometry necessary for mag- 
netic field generation. Astrophysical dynamos, however, 
are often generated by flows in simply-connected geome- 
tries which do not exhibit the scale separation between 
the large-scale magnetic field and the fluctuating part of 
the velocity field employed in the Riga and Karlsruhe ex- 
periments. This discrepancy has prompted the creation 
of several experiments to investigate the dynamo onset 
conditions in simply-connected flows with fully developed 
turbulence [11, 12, 13]. 

The threshold for the dynamo instability is not ex- 
pected to be the smooth transition from decaying to 
growing magnetic fields described by laminar kinematic 
theory [14]. Large-scale eddies can cause the instanta- 
neous flow to differ significantly from the mean flow. The 
growth rate of the magnetic field is highly sensitive to 
the flow geometry, and so the instantaneous flow may 
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FIG. 2: (a) PDF constructed from LDV measurements of 
V4, in the water model of the experiment (+) with Gaussian 
fit (dotted line). The difference between V4, and its mean is 
scaled to the standard deviation ay- (b) PDF constructed 
from Br measurements in the sodium experiment near the 
equator (+) and near the symmetry axis (o) with Gaussian 
fit (dotted line). The difference between Br and its mean is 
scaled to the standard deviation cts. 
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FIG. 3: Contours of Br{9,(f>) measured on the surface of the 
sphere. The applied field is subtracted from the measure- 
ments. This snapshot of the measured field corresponds to an 
induced dipole field transverse to the drive shaft axis. 



occasionally satisfy Rm > Rrriciit while the mean flow 
does not. The fastest-growing global magnetic eigen- 
mode would then fluctuate between growing and decay- 
ing states. The threshold of magnetic field growth there- 
fore has a range characterized by intermittent bursts of 
magnetic field growth. In this Letter the observation 
of an intermittently excited magnetic field in a simply- 
connected, turbulent flow of liquid sodium is reported. 
The structure of the excited field is consistent with the 
largest growing magnetic eigenmode predicted from a 
laminar kinematic model of the mean flow. 

The Madison Dynamo Experiment [Fig. 1] is a 
1 m diameter stainless steel sphere filled with liquid 
sodium [15]. Results presented in this Letter are from a 
turbulent two-vortex flow, similar to the flows described 
in [6], created by two counter- rotating impellers. The 
impellers are each driven by 75 kW motors and can 
achieve rotation rates up to 25 Hz, corresponding to 
■Rwtip ~ ^J■a<^LVtip = 150 based on the impeller tip speed. 

The mean flow is designed to generate growing mag- 
netic fields according to a laminar kinematic dynamo 
model [11]. According to the kinematic eigenvalue calcu- 
lations, for sufficiently large Rm the experimental flow 
is expected to excite a dipole magnetic field oriented 
transverse to the symmetry axis [Fig. 1]. An array of 
74 temperature-compensated Hall probes on the surface 
of the sphere provides measurements of the instantaneous 
multipole structure of the magnetic field induced by cur- 
rents in the liquid sodium. The Hall probes are capable 
of resolving changes in the magnetic field down to 0.3 G 
with a maximum range of ±170 G. A pair of magnetic 
field coils coaxial with the axis of rotation are used to 
apply a nearly-uniform 50 G axial field. This seed field 
brings the field induced by the flow above the noise level 
of the Hall probes. The applied field is sufficiently small 
that the strength of the Lorentz force is about 1% of fluid 
inertial forces. 



The velocity field of an identical-scale water model 
of the experiment is measured using Laser Doppler 
Velocimetry (LDV) [15]. The velocity measurements 
have Gaussian probability distribution functions (PDF) 
[Fig. 2(a)] as expected from a stationary turbulent flow 
according to the central limit theorem [16]. The mag- 
netic fluctuations measured near the axis of symme- 
try of the sodium experiment also have Gaussian dis- 
tributions [Fig. 2(b)]. In addition to these normally- 
distributed fluctuations, there are intermittent, large- 
amplitude magnetic bursts observed on probes near the 
equator of the experiment. The magnetic field dur- 
ing a burst has the spatial structure expected from the 
least-damped magnetic eigenmode from kinematic the- 
ory [Fig. 3]. The orientation of the transverse dipole is 
random for each burst so that the time-averaged induced 
field is axisymmetric. 
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FIG. 4: Time series of the energy in the transverse dipole 
field for an impeller rotation rate of 10 Hz. The diamonds 
mark the peak of a burst where the energy exceeds 50% of its 
maximum value. 

The bursts are ensemble averaged to determine typ- 
ical characteristics. A burst is defined to occur when 
the energy in the transverse dipole field exceeds a cer- 
tain threshold. For this analysis, the threshold is 50% 
of the maximum energy of the time series [Fig. 4]. This 
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TABLE I: The magnetic Reynolds number Rm based on the maximum speed in measured flows, duration of the measurement 
T, number of bursts A^;,, average burst rate ft, burst width ti,, growth rate At,, estimate of the overall fraction of time the flow 
is bursting rid, mean energy {E), and standard deviation of the energy as as a function of the rotation rate of the impellers 
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FIG. 5: The ensemble average of bursts from three time series. 
The averaged burst is used to calculate the growth rate and 
burst width in Tab. I. 



threshold is sufficiently small to capture a large number 
of bursts yet significantly larger than the mean energy 
(about two standard deviations above the mean energy 
for each time series). The bursts are averaged together 
and the growth rate is determined by an exponential 
curve fit [Fig. 5] . The results for various impeller rotation 
rates are reported in Tab. I. It should be noted that the 
strength of the fluctuations in the field is at most equal 
to the on-axis applied field strength of 50 G, hence the 
Lorentz force due to the fluctuations is weak compared 
to inertial forces. 

There are several possible mechanisms for the excita- 
tion of the transverse dipole field. Velocity field fiuctu- 
ations are large, with V/ (V) « 0.5 as determined from 
LDV measurements. These large fluctuations cause the 
peak flow speed to vary, which can be interpreted as vari- 
ation in Rm. Fluctuations at the largest scales can also 
cause changes in the shape of the flow leading to variation 
of i?TOcrit- The statistics of the small-scale fluctuations 
could change, also contributing to variation of Rrricnt- 
If the kinetic helicity of the small-scale eddies becomes 
sufficiently strong, the net current generation could give 
rise to the observed magnetic field bursts [17, 18]. Re- 
gardless of scale, subtle changes in the flow can adjust 
the instability threshold. 

To estimate the likelihood the flow is self-exciting, the 
kinematic growth rate from the velocity fleld measured 



FIG. 6: Kinematic growth rate versus Rm for the mean flow 
measured in the water experiment (solid) and an optimized 
flow (dashed) . The vertical lines identify Rmcrit for each case. 
The PDFs of Rm for flows with three different impeller rota- 
tion rates are shown to demonstrate the increasing overlap of 
the ranges of Rm and i?7Ticrit. 



in the water experiment is calculated [Fig. 6]. The solid 
line shows the results of the calculation using the mean 
flow, whereas the dashed line shows the growth rate cal- 
culated for an optimized flow geometry similar to and 
within the fluctuation levels of the measured flow [14]. 
The region between the growth rates of the mean and op- 
timized flows indicates the range of possible eigenmode 
growth rates for an instantaneous realization of the flow 
and the resulting variation in Rnicrit- 

The velocity fluctuations also contribute to variations 
in the instantaneous maximum speed of the flow, thereby 
creating a range of Rm. Probability distributions of 
Rm constructed from the measured velocity fluctuations 
for three different impeller rotation rates are plotted in 
Fig. 6. It is expected that a greater overlap of the PDF of 
Rm with the range of Rmcrit will result in magnetic field 
excitations with greater frequency and strength. The du- 
ration of each excitation is expected to decrease since 
the correlation time of the velocity fluctuations scales as 
Tc ^ i/Ve where i is the eddy scale length and Vg is 
the characteristic speed of the eddy. For example, LDV 
measurements from the water model give Tc = 80 ± 20 
ms for the /tip = 16.7 Hz fiow, consistent with eddies of 
size i = 0.25 m and speed Vi — 3m/s. The proportion 
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FIG. 7: The PDF of the energy in the transverse dipole field 
for several impeller rotation rates. The thick line is an exam- 
ple Poisson distribution. The thin line represents the energy 
distribution if the magnetic fiuctuations were Gaussian. 



of time that the magnetic field is bursting is estimated 
to be Ud = fbTb, where fb is the average frequency of the 
bursts and is the width of the conditionally-averaged 
burst at half-maximum. The data in Tab. I show that 
the proportion of time the flow is bursting stays relatively 
constant between 5-8%. 

Table I reveals that the standard deviation of the en- 
ergy in the intermittent transverse dipole field is approxi- 
mately equal to its mean value, a characteristic of a Pois- 
son probability distribution [19]. Assuming each excita- 
tion can be treated as a rare random event, the prob- 
ability distribution of the magnetic field energy can be 
determined heuristically. The probability of measuring n 
bursts in time t is given by P{t) ~ {fbt)'^e~^**/n\ where 
/b is the average rate of bursts. The average growth of the 
magnetic field over time t during a burst is AS = Bge^*, 
where Bq is the average strength of the initial seed field. 
The resulting gain in the magnetic field energy per unit 
volume is A£' = AB^/2fiQ = (S^/2^o) exp(2At) and so 
t = \og{AE/Eo)/2\ where Ea = B^/2^o- Substituting 
the time in terms of AE into the Poisson distribution 
yields a log-Poisson distribution for the probability den- 
sity of AE: 



P{AE) 



1 



2A V ^0 



-(/i,/2A)ln(A£;/£;o) 



(2) 



The probability distributions of the transverse dipole en- 
ergy are shown in Fig. 7. The distributions with large 
numbers of bursts tend to have significantly more high 
energy fluctuations than is expected from Gaussian fluc- 
tuations. The overall invariance of the distributions as 
the impeller rotation rate is increased demonstrates that 
the increased frequency of bursts is offset by their short- 
ened duration. 

The results presented demonstrate how turbulence in 
a simply-connected geometry changes the onset condi- 
tions of the dynamo. Rather than the smooth transi- 
tion from damped to growing fields predicted by either 
kinematic or mean field dynamo theory, the transition is 



characterized by intermittent magnetic field bursts which 
may be relevant to some dynamo models [20]. Although 
sustained growth is not yet observed, the transient exci- 
tation demonstrates the intermittent characteristics of a 
turbulent dynamo. 
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